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We present results of theoretical study and numerical calculation of the dynamics of molecular 
liquids based on combination of the memory equation formalism and the reference interaction site 
model - RISM. Memory equations for the site-site intermediate scattering functions are studied in 
the mode-coupling approximation for the first order memory kernels, while equilibrium properties 
such as site-site static structure factors are deduced from RISM. The results include the temperatu- 
re-density(pressure) dependence of translational diffusion coefficients D and orientational relaxation 
times t for acetonitrile in water, methanol in water and methanol in acetonitrile, all in the limit of 
infinite dilution. Calculations are performed over the range of temperatures and densities employing 
the SPC/E model for water and optimized site-site potentials for acetonitrile and methanol. The 
theory is able to reproduce qualitatively all main features of temperature and density dependences 
of D and r observed in real and computer experiments. In particular, anomalous behavior, i.e. the 
increase in mobility with density, is observed for D and r of methanol in water, while acetonitrile 
in water and methanol in acetonitrile do not show deviations from the ordinary behavior. The 
variety exhibited by the different solute-solvent systems in the density dependence of the mobility 
is interpreted in terms of the two competing origins of friction, which interplay with each other as 
density increases: the collisional and dielectric frictions which, respectively, increase and decrease 
with increasing density. 
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I. INTRODUCTION 

The mobility of either polar molecules or small po- 
lyatomic ions in dense dipolar liquids exhibits some 
anomalies that have attracted considerable interest of 
scientista 1 iVi 4 i5i6,7,8,V°i 11 i 12 , 1 3, 14 i^ In this paper we 
pay attention to the study of the so-called "anomalous 
pressure behavior" of polar binary molecular liquids by 
considering such popular models as methanol in water, 
acetonitrile in water and methanol in acetonitrile. In 
particular, we are interested in the dcnsity(pressure) de- 
pendence of translational diffusion coefficients and reori- 
entation relaxation times of solutes. To begin with, we 
would like to say first few words about anomalies as they 
are. 

The very primitive (classical, not quantum) picture of 
a molecular liquid could be the set of finite size struc- 
tural particles - molecules - consisting of interacting 
sites, which are connected by chemical bonds. Positions 
of molecules as well as their orientations are not fixed in 
general and may vary randomly in time. The change may 
be achieved either by translational or rotational motion. 



Intuitively, it is then clear to understand that applying 
the external pressure for such system results in making 
the intermolecular distances shorter, which as the con- 
sequence will limit the molecular mobility. In this case 
one may expect a monotonous decrease with pressure for 
such quantities as translational or rotational diffusion co- 
efficients. Limitation of the molecular mobility can also 
be discussed in terms of some other dynamical charac- 
teristics such as viscosity, dielectric relaxation, etc. This 
point of view has been confirmed experimentally for nu- 
merous liquids and solutions^. At the same time, some 
molecular liquids do not obey strictly this abstraction 
and exhibit behavior which is different from the intuitive 
picture. Typical representatives of such molecular liquids 
are water and some aqueous solutionsiiifiiiS whose mo- 
bility with pressure can be enhanced. The enhancement 
of the molecular (or ionic) mobility under pressure then 
has been identified as the anomalous one. 

Historically, first attempts to explain dynamics of 
molecules in solutions were based on the use of the so- 
called continuum model for solventa 20 i 21 i 22 i 23 i 24 i 25 . In 
this model the surrounding medium is regarded as a di- 
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electric continuum extending up to the surface of the par- 
ticle. The particle itself is treated as a rigid sphere of 
finite size with some charge for ions or dipole moment 
for molecules. Dynamic behavior of the media are de- 
scribed by the Debye equation in terms of the frequency 
dependent dielectric susceptibility and hydrodynamic ef- 
fects are entirely neglected. The most important contri- 
bution of the theory was that it succeeded in extending 
the phcnomenological approaches to consider the friction 
on a rotating dipole. When a dipole is surrounded by 
a dielectric continuum, it polarizes the medium. If the 
dipole rotates, the polarization should relax so that it 
adjusts itself to the electric field produced by the new 
orientation of the dipole. The relaxation process is asso- 
ciated with the energy dissipation which is identified as 
the dielectric friction. Later the theory was improved by 
obtaining velocity of fluid from linearized Navier-Stokes 
equation for a viscous incompressible fluid with bound- 
ary conditions of either perfect sticking or slipping at the 
surface of the particle2£i2L2&. Subsequently, there were 
many another suggestions on improvements of the theory, 
e.g. the enhanced "continuum mechanical framework" 
description has been given in series of works22i22i2i, where 
dissipation function is constructed by adding the electri- 
cal dissipation to the ordinary hydrodynamic one, etc. 
However, contrary to what was expected, the contin- 
uum model turned out to be unable to explain properly 
numerous contradictions of phcnomenological laws with 
experimental data. One example is the failure of the 
Stokes-Einstein relation in explaining the dependence of 
mobility of ions in polar liquids on their size. In or- 
der to compromise the trouble it has been suggested to 
use the so-called "effective" radii of particles, or rewrite 
the relation so that it appears in the fractional exponent 
formS. Such kind of revisions had some effects, but there 
were also other difficulties which could not be resolved 
in principle. For instance, dielectric continuum theories 
cannot explain observed difference in mobility of anions 
and cations because the interaction between an ion and 
dielectric media does not depend on the sign of the ion 
charge. That is enough to realize that the insolvency of 
a conventional continuum model description can be over- 
come only by a microscopic formulation of the problem. 

The microscopic, or statistical-mechanical, theory of 
molecular liquids is now on its rise, and the consider- 
able progress has been achieved during last years both 
in the description of equilibrium and non-equilibrium 
properties^. It became possible due to the aptly 
used reference interaction site model - RISM 34,35,36 
and its modifications, or in the case of dynamics - 
due to appropriate unification of RISM and general- 
ized Langevin equation (GLE) which is obtained in the 
Mori formalism 37 - 38 .. The choice in favor of generalized 
Langevin equation is not accidental. For example, an al- 
ternative way of describing non-equilibrium processes can 
be based on the use of various types of kinetic equations 
which are derived from the BBGKY hierarchy 39 . How- 
ever, when it is applied to systems with high density, 



there is an ambiguity in the construction of the collision 
integral. Also the linearization procedure for the colli- 
sion integral cannot be justified by the same arguments 
which are used in the case of systems with low density 4 ^. 
In Mori formalism such questions do not appear. GLE 
describes the time evolution of several dynamic variables 
in concern as a function of the representative point in 
the phase space. All other variables which are not un- 
der explicit concern are projected onto the dynamic vari- 
ables of interest with the help of a projection operator. 
The projection leads to an equation which looks similar 
to the phenomenological Langevin equation containing 
the factional force proportional to the rate of change 
of the dynamic variables, and the random force, which 
are related to each other by the fluctuation-dissipation 
theorem. If one chooses as the dynamic variables the 
density and the conjugated current of sites or atoms of 
molecules in liquid, the theory gives a memory equa- 
tion for the dynamic structure factor of atoms, which 
describes the time evolution of the site-site density pair 
correlation functions 4 - 1 . Initial values, i.e. site-site static 
structure factors and direct correlation function which 
are necessary to solve this integro-differential equation, 
are obtained from the RISM or one of its generalizations. 
A crucial development of the theory is rather conceptual, 
not mathematical, in the sense that it has provided a new 
concept to view dynamics of a molecule in solution, which 
is quite different from the model traditionally exploited 
in the field. The new model regards the molecular motion 
in liquid as a correlated translational motion of atoms: 
if two atoms in a diatomic molecule are moving in the 
same direction, then the molecule as a whole is making 
translational motion, while the molecule should be rotat- 
ing if its atoms are moving in opposite directions. This 
view differs from the traditionally developed rotational- 
translational model which is based on the use of angu- 
lar coordinates^2iM44iiSii&. Another important issue of 
RISM-based theories, which should be mentioned here, 
is that they pioneered in explaining the observed differ- 
ence in the behavior of anions and cations in solution. 
Both equilibrium^ 7 " and non-equilibrium 4 ^ systems were 
studied, and now we know that these differences should 
be attributed to dissimilarities of the microscopic struc- 
ture of solvent surrounding the ion. 

The new theory of dynamic processes in molecular liq- 
uids has been under development in last decade in works 
by Hirata and Chong 4 ^*^ and get its extension in studies 
by Hirata, Yamaguchi and othera 51 i 52 i 53 i 54 i 55 i 56 i 57 i 58 i 59 i 60 . 
Particularly, theoretical description of pressure depen- 
dence of the molecular motion of liquid water and trans- 
port properties of ionic liquids were subjects for inves- 
tigation in recent papers by Yamaguchi et alm^&. In 
present manuscript we continue this line and report our 
latest study which is the pressure and temperature de- 
pendence of dynamics of polar molecules in solution. The 
manuscript is organized as follows. In subsection III Al 
we write down equations of motion, which are memory 
equations, for site-site intermediate scattering functions, 
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and give relations for translational diffusion coefficient 
and reorientation correlation time in terms of these func- 
tions. In subsection III Bl we deal with equilibrium prop- 
erties of the system which we obtain using DRISM with 
hypernetted chain (HNC) type of the closure. In order 
to solve a memory equation one needs to set the model 
for memory kernels. This question is discussed in some 
detail in subsection III CI More attention is paid to the 
so-called mode-coupling approximation, which is used to 
study dynamical behavior of the system in the long-time 
scale. We apply our formalism to investigate some realis- 
tic systems such as acetonitrile in water, methanol in wa- 
ter and methanol in acetonitrile. Set up for these models 
and details of numerical procedures in the computational 
part of this work are briefly explained in section ITTT1 Ob- 
tained results and remedies are discussed in section Hvl 
and conclusions are given in section Ivl 



II. THEORY 



Memory equations for site-site intermediate scattering 
functions 



Time-correlation functions are convenient tools for de- 
scribing properties of many-body systems and can direct- 
ly be associated with experimental observables. A gene- 
ral theoretical scheme for calculation of time-correlation 
functions is provided by the Mori approach 3 ^^, where 
they are phrased in terms of memory function equations. 
Hereafter we will follow the formalism by Hirata and 
Chong 50 , which is the judicious unification of the the- 
ory of dynamical processes in simple liquids based on the 
Mori approach 37 - 38 , and reference interaction site model 
for molecular liquids 3 ^ 3 ^ 3 ^. In this formalism "slow" 
variables of the system are partial number densities and 
longitudinal current densities. Then, of practical inter- 
est are elements of the matrix of the site-site intermediate 
scattering functions F(fc; t) and the matrix of their self- 
parts F s (k;t) defined by 



F^{k;t) = i(p«'*(k;0)pT(k;t)), (2.1a) 
Fr(k;t) = i(^'*(k;0)p^(k;t)) s , (2.1b) 



where N is the total number of particles, p a (k;t) is the 
site a number density in reciprocal space, k is the wave- 
vector, k = |k|, t is time, * means complex conjugation, 
angular brackets denote appropriate statistical average 
(e.g., canonical Gibbs ensemble average), and suffix "s" 
stands for "self" and means the correlations between two 
sites in a same molecule. In the present paper, we con- 
sider the limit of infinite dilution, the limit that the solute 
concentration vanishes. One has two types of memory 
equations for the solvent subsystem and one for the so- 
lute subsystem (indicated by the superscript "u" ) which 



arc: 



F(k;t) 



F-(M) 



F.(M) 



-<w£)F(fc;t)- / dTK(fc;r)F(M-r), 
Jo 



(2.2a) 



ul) t F B {k;t) - drK s (k;T)F s (k;t-T) 



(2.2b) 

-<u, 2 fe ) s u F s u (fc; t) - f At K s u (fc; r)F s U (fc; * - r). 
Jo 

(2.2c) 



In these equations, dot over the quantity means its time 
derivative, and the memory function matrices, denoted as 
K(k;t), K s (k;t) and K"(fc;t), describe the friction on the 
motion of interaction sites and are subject for definition 
later. Quantities (w?) a and (w 2 ,)" are normalized 

second order frequency matrices given by 



,2\u 



fc 2 J L (fc)S^ 1 (fc), (2.3a) 
k 2 Uk)S-\k), (2.3b) 
k 2 mS^-\k), (2.3c) 



where S(fc) = F(k;t = 0), S s (k) = F s (k;t = 0) and 
Sg (k) = Fg (k; t = 0) are matrices of static site-site struc- 
ture factors and their self parts, respectively, while J L (fc) 
is the matrix of static longitudinal site-current correla- 
tion functions, i.e. 



N ^ 



ik- (r 



(2.4) 



where subscripts i,j refer to molecules, rf = rf (0) is the 
initial position of site a, and vf z = vf z (0) is longitudinal 
component of the initial velocity of site a. The analytical 
expression for arbitrary shape of the molecule has been 
given recently by Yamaguchi et al^> and reads 



LJ L (fc)j^: 

k B T 



knT 



(#*7)2 



\jo{kl^)+h{kl a ^)] 
x[«5r a ] T - [Tr l _1 l - I -1 ] • <5r 7 
j 2 (k£ a ~<)[Sr a x 5r~<} T ■ I 1 • [5r a x <5r 7 ], 

(2.5) 



where fee is the Boltzmann constant, T is thermodynamic 
temperature, M is the total mass of the molecule, l ai is 
the distance between sites, 5r a is the vector pointed from 
the center of mass to the site a, I is tensor of inertia 
moments of the molecule, 1 is the diagonal unit matrix; 
finally, jo and ji are spherical Bessel functions of the first 
kind. Definition for J"(fc) is similar to the one given by 
equation l|2.4|l with the difference that summation runs 
over the solute molecule only. 
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Expressions for static site-site structure factors and 
memory kernels are discussed separately in sections IIIII 
and II VI respectively. Here we show how quantities of our 
interest in this paper, which are solute's translational dif- 
fusion coefficient D and rank-1 reorientation correlation 
time tri, are related with F^(k;t). The diffusion coef- 
ficient can be expressed in terms of the site-site veloc- 
ity autocorrelation function, while the relaxation time is 
integrated rank-1 reorientation autocorrelation function. 
Both of them are described in terms of the F"(k; t). 

Following the derivation procedure presented in Ref. 
IrUl the translational diffusion coefficient D is obtained as 



In this work we restrict our consideration to reorienta- 
tion relaxation processes of the rank-1 only and hereafter 
will not use the subscript Rl in notations. For the (rank- 
1) reorientation relaxation time we use definition^ 



dtC^t). 



(2.10) 



D 



dt Z a ~>(t) 



1 r°° i 



= — lim 



^l^[F s (fc;r)j»\ 



(2.6) 



where Z ai (t) is the site-site velocity autocorrelation 
function with sites a and 7 belonging to the same 
molecule. 

The rank-1 reorientation autocorrelation function 
Cfj,(t) is defined by 



(/*«(Q)/ii(*)) 
£,-<N 2 > ' 



(2.7) 



where is a vector fixed on the molecule i. In our case 
it is the dipole moment and therefore can be described 
by the linear combination of site coordinates as 



B. Equilibrium properties of molecular liquids via 
RISM/DRISM 

Initial values of intermediate scattering functions, that 
we need to solve memory equations 1)2. 2p . can be ob- 
tained using the RISM theor}r 33 i 34 i 35 i 36 . It predicts static 
structure of molecular fluids via the calculation of site- 
site pair correlation functions. This method has been 
extensively used and proved to be the powerful tool in 
the microscopic description of equilibrium quantities of 
the system 3 ^. In the present article, we employ the 
DRISM version of the theory 64,6 - instead of the standard 
RISM/HNC theory. The DRISM theory incorporates a 
bridge function which is determined to reproduce the di- 
electric constant of solvent in phenomenological fashion. 
The use of the DRISM theory is known to improve the di- 
electric loss spectrum of solvent quantitatively, although 
it does not change the essential conclusion of the study 
with respect to the pressure dependence of the dynamics 
of solute in solvent™. 



/**(*) 



z a v a (t) 



(2.8) 



with z a being site partial charges. Since the molecule is 
considered to be the neutral one, partial charges satisfy 
the condition of electro-neutrality ^ Q z a = 0, and make 
definition of /x,(i) invariant under the transformation of 
translation of the laboratory reference frame. Putting 
equation (|2.8|) into equation (|2.7|) and using properties of 
time-correlation functions 62 one arrives at 



£i£« 7 *a*7 (r?(0)r7(*)) 

£,<N 2 ) 

Nj: aj z a z 7 Z a ~f(t) 

£;(N 2 > ' 



(2.9a) 
(2.9b) 



Hence, the time development of both Z ai (t) and C^t) 
is governed by the memory equation for the self-part of 
the site-site intermediate scattering function F"(fc;£). 



C. Mode-Coupling approximation for memory kernels 

It makes sense to talk about translational diffusion co- 
efficient and reorientation relaxation time only when dy- 
namics of the system is considered in the long-time limit, 
i.e. when the time scale is large enough for fast relax- 
ation processes to be completed. Memory kernel of the 
memory equation in that case is usually constructed by 
the so-called mode-coupling approximation 63 . In works 
by Chong et alm^L the conventional mode-coupling the- 
ory has been extended to the case of molecular liqu- 
ids based on the interaction-site model. In particular, 
mode-coupling expressions for memory function matri- 
ces, which we will denote as K MC (k,t) and K S>MC (&, t), 
were obtained as bilinear combinations of site-site inter- 
mediate scattering functions or their self parts as 
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LJ^(*)-Kmo(M).T = j£yt y dq {<gL£(9) * F(«.*) ' c( 9 )J«TF^(|k- q|,t) 

- q z (k - q z )lc(q) ■ F(<j, t)J Q7 [F(|k - q|,t) • c(|k - q|)j^}, (2.11a) 
[J-W-KsMk^r = 7^)3 /dqg 2 2 L£(<z)-F(g,t).£(g)j^Fr(|k-q|,i), (2.11b) 

here for simplicity the wave-vector k has been taken to be parallel to z-axis. It has been shown, however, that 
the proposed expressions for memory functions underestimate friction in orientational motions^. According to the 
recipe by Yamaguchi and HirataS*, memory functions for the self-part should be given by the linear combination of 
corresponding mode-coupling memory functions as 

ik„(m) ■ uk)r = £ ^ 1 zz m2 zz m A^' (i? - T ^ Tl+T ' ) ) ia-'w ■ iwmxt, (2.12) 

mi,2,3 = {a:,j/,z} 

(i) 

where u\m stands for the unitary matrix that describes rotation between molecular and space-fixed reference frames 
of the molecule i, while Z^[ m2 denotes orientation-dependent site-site velocity correlation matrix in the molecular 
coordinate system and is given by 

7 rr\ 

Z Zm 2 = -jfSmum, + x 5r a ] T ■ I -1 ■ [e„ l2 x 6tT]. (2.13) 

Here e TOl is the unit vector in the m\ direction. The collective part of the memory function (neglecting the orientational 
correlation between different molecules) is obtained as^ 

K(fc,t) = K MC (fc,<) + K,(fe,t) - K s , M0 (fc,t). (2.14) 
In the case of solute corresponding quantities are given by relations 

LJr^K s u MC (M)J Q7 = J^fj dqgfL£ UV (9)-F v ( g> t)-c vu ( g )j^Fr 7 (|k-qM), (2.15) 

"H, 2l3 ={x,3/,z} 

Here, superscript "v" is used to indicate the solvent subsystem, and intermediate scattering functions F v (fc,i) are 
those obtained in the process of solving of memory equations l|2.2a[l and (|2.2b|l with mode-coupling expressions for 
memory kernels deduced from equations (|2.14(1 and l)2.12[l . respectively. 



III. SETUP OF MODELS AND NUMERICAL 
PROCEDURES 

We performed explicit calculations for several popu- 
lar systems, namely, acetonitrile (CH3CN) in water, 
methanol (CH3OH) in water and methanol in acetoni- 
trile, all in the case of infinite dilution. As for the 
structure and the intermolecular potential of water we 
employed a model of the extended simple point charge 
(SPC/E)S£. We also put the Lennard- Jones (LJ) core on 
the hydrogen atoms in order to avoid the undesired diver- 
gence of the solution of the RISM integral equation. The 
LJ parameters of the hydrogen atom, the depth of the 
well and the diameter, were chosen to be 0.046 kcal/mol 
and 0.7 A, respectively. 

In acetonitrile and methanol the methyl group was con- 



r 



sidered to be a single interaction site (Me) located on the 
methyl carbon atom. Consequently, all chemical com- 
pounds consist of three sites which interact through the 
pair potential§2i2!i 



_ 



— 4£an 



' 07 



; Q7 



+ 



■^0^7 



(3.1) 



i.e., LJ plus Coulomb. Here r*? 



r 7 |; parameters 



ai and a ai are LJ well-depths and LJ diameters defined 



e 

as e 



c*7 



^/e Q e 7 and a aj = (a a + <7 7 )/2, respectively, 
with a a being the LJ diameter of a single site. Point 
charges for acetonitrile were chosen to reproduce elec- 
trostatic potential obtained in ab initio calculations^. 
Numerical values of masses of sites, parameters of the 
site-site interaction potential (13.111 and cartesian coordi- 
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nates of sites are specified in Table [H Information about 
bond length can be deduced from cartesian coordinates 
of sites (x, y, z). 



TABLE I: Masses of sites m a , parameters of the site-site in- 
teraction potential 13.111 and cartesian coordinates of sites for 
watei^, acetonitrile^ and methanol. All values of elj are 
multiplied by the factor 10 14 , a.u. stands for "atomic units". 



stant for acetonitrile at different temperatures are indi- 
cated in Table |ni| 





m a 




Co 




X 


y 


z 




(a.u.) 


(a.u.) 


(A) 


(erg/molec) 


(A) 


(A) 


(A) 


o 


16.0 


-0.8476 


3.16 


1.084 








-0.0646 


H 


1.008 


0.4238 


0.7 


0.3196 





0.8165 


0.5127 


H 


1.008 


0.4238 


0.7 


0.3196 





-0.8165 


0.5127 


Me 


15.024 


0.269 


3.6 


2.64 








1.46 


C 


12.0 


0.129 


3.4 


0.6878 











N 


14.0 


-0.398 


3.3 


0.6878 








-1.17 


Me 


15.024 


0.265 


3.74 


1.4525 


-1.4246 








O 


16.0 


-0.7 


3.03 


1.1943 











H 


1.008 


0.435 


1.0 


0.3196 


0.3004 


0.8961 






In calculations for acetonitrile or methanol in water 
the temperature of the system was varied from 258.15 to 
373.15 K, and the density of water from 0.9 to 1.2 g/cm 3 ; 
for the case of methanol in acetonitrile the temperature 
of the system was varied from 293.15 to 323.15 K, and 
the density of acetonitrile from 0.6726 to 0.815 g/cm 3 . 
Connection of the water parameters with thermodynamic 
pressure is shown in Table [H] (except for the metastable 
regions where we do not have reliable data). 



TABLE II: Density-pressure correspondence for wates^ given 



for temperatures Ti 
373.15 K. 



273.15 K, T 2 = 298.15 K and T 3 



Density 


Number density 


Pressure at Ti 


Pressure at T2 


Pressure at T3 


(g/cm 3 ) 


(A" 3 ) 


(MPa) 


(MPa) 


(MPa) 


0.900 


0.03008 








1.000 


0.03334 


0.4085 


6.6914 


100.6450 


1.025 


0.03426 


52.6119 


66.3549 


173.0629 


1.050 


0.03510 


111.7871 


133.7116 


255.2558 


1.075 


0.03593 


179.4861 


209.9492 


347.8493 


1.100 


0.03676 


257.2009 


296.1955 


451.5134 


1.125 


0.03760 


346.1457 


393.4575 


566.9920 


1.200 


0.04011 


689.0266 


760.7554 


993.3837 



TABLE III: Density and dielectric constant for acetonitrile 
as functions of temperature: experimental data used in our 
computation. Temperature is in K, and density is in g/cm 3 . 

T 293.13 295.05 298.15 303.15 308.15 313.15 323.15 
~~p .782 a — .7762 b .7712 a .7652 b .7603 a .7492 a 
e 38. 8 C 37.5 d 36.69 c 35.93 c 33.5 f 

a Rcf. [zl b Ref. |2l c Ref. d Ref. [ftj e Ref. \zH f Extrapolated. 

From the memory equation / mode-coupling theory 
and the DRISM/HNC integral equation theory, the diffu- 
sion coefficients and the reorientation relaxation times of 
solute molecules in solution can be obtained based solely 
on the information about molecular shapes, inertia pa- 
rameters, intermolecular interaction potentials, temper- 
ature and density. First, we calculate the site-site static 
structure factor by solving the DRISM equation using 
the intermolecular interaction, molecular shape, tempe- 
rature and density. In order to improve the convergence 
of the DRISM calculation, we use the method of the 
modified direct inversion in an iterative space (MDIIS) 
proposed by Kovalenko et ali& From the static site-site 
structure factor, we calculate the site-site intermediate 
scattering function using the site-site memory equation 
with the mode-coupling approximation for the memory 
kernels. The memory equation is time-integrated nume- 
rically. Time-development of correlation functions in the 
k — > limit is treated separately by the analytical limi- 
ting procedure of theoretical expressions. In the numeri- 
cal procedure, the reciprocal space is linearly discretized 
as k = (n+ |)Afc, where n is an integer from to Nk — 1. 
Values of Afc and Nk are 0.061 A -1 and 2 9 = 512, respec- 
tively. The choice for Nk as the power of two is coming 
as the requirement of the subroutine for the fast Fourier 
transform, which is used in DRISM / MDIIS. The diffu- 
sion coefficient D is calculated from the asymptotic slope 
of the time dependence of the mean square displacement, 
and the orientational relaxation time r is deduced from 
the rotational autocorrelation function. 



IV. RESULTS AND DISCUSSIONS 



Temperature/density dependent dielectric constant e 
for water used in numerical calculations has been eval- 
uated as a physical solution of an empirical nonlinear 
equation presented in Ref. 



e — 
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9.32 ■ 10 4 (1 + 



153 ' 
i-T 1 / 4 / 



(1-3/ufT 



(3.2) 

where v is a molar volume in units of cm 3 /mol, and T 
is thermodynamic temperature in K. This equation was 
also used in such temperature/density points where no 
experimental values exist. Density and dielectric con- 



A. Equilibrium properties: static structure 

Since the discovery of anomalous behavior of the mo- 
lecular mobility with pressure it has been attributed to 
the hydrogen bonding properties of the system. The idea 
is that the hydrogen bonds or hydrogen bonding net- 
work are distorted upon compression, so that the dis- 
solved molecule is actually having fewer hydrogen bonds 
at higher pressure, thus making it easier to move and ro- 
tate. The extensive test of this idea has been the objec- 
tive of many researches. In particular, the time evolution 
of the formation and rupture of hydrogen bonds turns 
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2.4 



FIG. 1: Site-site radial distribution functions of neat water at 
T = 273.15 K and set of densities from 0.9 to 1.125 g/cm 3 , ob- 
tained by the DRISM/HNC integral equation theory. Arrows 
show directions of alternations due to an increase in pressure. 



out to be a popular quantity to be examined in studies 
based on the computer simulation (see, e.g., Ref. l79l and 
references therein). In order to circumvent the problem, 
a number of alternative definitions of hydrogen bonds 
was suggested on top of the traditional one based either 
on geometrioStt and/or energetio^i criteria, though none 
of them is conclusive. The definition of hydrogen bond 
in the theoretical study based on the integral equation 
is somewhat ambiguous as well. However, it is a quite 
general observation for such liquids like water and alco- 
hol that a pair of atoms having small size and/or large 
charges with opposite sign makes a well-defined first peak 
in the pair correlation function: for example, hydrogen 
and oxygen in water and alcohol. Especially, in order 
for making a hydrogen bond, one of the pair of atoms 
should be as small as a hydrogen atom. Since those li- 
quids are known to make a hydrogen bond between the 
pair of atoms, the well-defined peak has been assigned to 
a "hydrogen bond" . The hydrogen bond between a pair 
of atoms induces a "liquid structure" which is specific 
to molecules consisting the fluid, which in turn produces 
peculiar peaks in the pair correlation function at parti- 
cular distances. For example, in the case of water, the 
hydrogen bond peak appears at ~ 1.8 A in the O-H cor- 
relation function, which induces a new peak in the 0-0 
correlation function around ~ 4.5 A. The 0-0 distance is 



O 1.2 
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FIG. 2: Water-methanol site-site radial distribution functions 
for sites, that may be related to making hydrogen bonds, at 
T = 273.15 K and set of densities from 0.9 to 1.125 g/cm 3 
for water, obtained by the DRISM/HNC integral equation 
theory. In the notations used, first site always belongs to 
water, and second site always belongs to methanol. Arrows 
show directions of alternations due to an increase in pressure. 



peculiar to the ice-like tetrahedral coordination. There- 
fore, the peak at ~ 4.5 A is regarded as a "finger print" 
of the hydrogen bond network existing in liquid water—. 
Alcohol also makes a hydrogen bond between a pair of 
the molecules, but does not make a three dimensional 
network, which is distinct from water. Instead, it makes 
zigzag chains of hydrogen bonds. It is because one of the 
four hydrogen-bonding sites is replaced by an alkyl group 
which is too large to form the bond. In the case of neat 
acetonitrile, there is no indication for making hydrogen 
bonds between a pair of the molecules. It is understood 
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FIG. 3: Water-acetonitrile site-site radial distribution func- 
tions at T = 273.15 K and set of densities from 0.9 to 1.125 
g/cm 3 for water, obtained by the DRISM/HNC integral equa- 
tion theory. In the notations used, first site always belongs to 
water, and second site always belongs to acetonitrile. Arrows 
show directions of alternations due to an increase in pressure. 



FIG. 4: Acetonitrile-methanol site-site radial distribution 
functions at T = 293.15 K and set of densities from 0.6726 to 
0.815 g/cm 3 for acetonitrile, obtained by the DRISM/HNC 
theory. In the notations used, first site always belongs to ace- 
tonitrile, and second site always belongs to methanol. Arrows 
show directions of alternations due to an increase in pressure. 



from its molecular constituents, i.e. it does not have a 
small atom as a possible hydrogen-bonding site. There- 
fore, the structure of the liquid or the pair correlation 
function is determined essentially by the repulsive inter- 
actions and the molecular geometry. 

The pressure or density dependence of the liquid struc- 
ture is determined essentially by the interplay between 
the repulsive interaction and the hydrogen bonding, 
which interfere each other through the intramolecular 
constraints or molecular geometry. Increasing pressure 
will make repulsive interactions more dominant, which 
will weaken the hydrogen bond by distorting the favo- 
rable arrangement between the two molecules. The pair 
correlation functions on pressure or density should re- 
flect the interplay mentioned above. Increasing pressure 
or density will result in reduction in all the indications of 
hydrogen bonding and of the consequential liquid struc- 
ture, and will make the steric or packing effect more do- 
minant. 

Fig. Q] shows the site-site radial distribution functions 
(RDF's) for neat water. The sharp peak at ~ 1.8 A in 
the O-H RDF, Fig. [TJd), is that of the hydrogen bond 
without any question. A weak but conspicuous second 
peak at ~ 4.5 A in the 0-0 RDF, Fig. [IJl), is indica- 
tive of the tetrahedral configuration characteristic in the 



hydrogen-bond network in water. All the changes in the 
RDF caused by increasing density coincide with the intu- 
itive picture emerged from the disruption or distortion of 
the hydrogen bond and its network: the reduction of the 
hydrogen bond peak in the O-H RDF, the reduction of 
the second peak in the 0-0 RDF. The increasing height 
in the first and third peaks in the 0-0 RDF implies the 
enhanced tendency of the steric and packing effect due 
to the increased density. 

Fig. |21 shows the site-site RDF for methanol in water. 
There are two possible hydrogen bond sites in methanol, 
O and H, which make hydrogen bonds, respectively with 
H and O of water molecules. The 0(water)-H (methanol) 
RDF, Fig. has the hydrogen bond peak at the 

position where that in neat water appears. The peak 
height and resolution are less compared to the case in 
neat water, implying the hydrogen bond is somewhat 
weaker. The density dependence of the peak shows si- 
milar behavior with that in neat water, indicating that 
the hydrogen bond is weakened due to the steric dis- 
tortion caused by pressure. The 0(methanol)-H(water) 
RDF, Fig. 13;), has the hydrogen bond peak at ~ 1.5 
A, which is much shorter than that in the O(water)- 
H(methanol) RDF. The peak is even more sharply de- 
fined than that in neat water, implying that the hydro- 
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gen bond is stronger. The density dependence of the 
peak shows similar tendency with that in neat water, 
but decreases more sharply with increasing density. The 
0(water)-0 (methanol) RDF, Fig. QJi), maintains the es- 
sential structure which is seen in that of neat water. The 
important point is that the second peak at ~ 4.5 A has 
not disappeared, but is even enhanced, which indicates 
the tctrahcdral coordination peculiar to the hydrogen 
bond network in water keeps its integrity. Indeed, in our 
earlier work concerning the structure of the tert-butyl 
alcohol-water mixture, we have clarified that a butanol 
molecule is incorporated in the hydrogen-bond network of 
water in its infinite dilution by forming hydrogen bonds 
with surrounding solvent molecules^. 

The site-site RDF of the acetonitrile in water is shown 
in Fig. [3]along with its density dependence. According to 
the criteria for hydrogen bonding ability of atoms stated 
above, the acetonitrile molecule has one possible hydro- 
gen bond site, nitrogen N, which can act only as a pro- 
ton acceptor. As is anticipated from the argument, the 
N(acetonitrile)-H(water) RDF has a conspicuous hydro- 
gen bond peak at about ~ 1.8 A, which again decreases 
with increasing density, Fig. Ob)- On the other hand, 
the N(acetonitrile)-0 (water) RDF shows less conspicu- 
ous second peak at 4.5 A, which implies that the hydro- 
gen bond network of water is somewhat disturbed around 
the solute compared to the neat solvent, Fig. 

In Fig. 01 shown is the site-site RDF for methanol in 
acetonitrile at infinite dilution. In this solvent 
(acetonitrile) molecule can make a hydrogen bond with 
the solute molecule only at the N-site by accepting the 
proton from methanol. The first peak in Fig. 0Jd) is 
indicative of the hydrogen bond. The solvent molecules 
do not make any hydrogen bond among themselves. The 
N-site never forms a hydrogen bond with the oxygen site 
of methanol. It is not only because the both sites have 
negative charges, but also because the neither of them 
is as small as hydrogen. Therefore, the N(acetonitrile)- 
O (methanol) RDF is largely determined by the repulsive 
core and geometric constraints. The system is a good 
example of the case in which constituent molecules can 
make only one hydrogen bond, only between solute and 
solvent. 



B. Dynamical properties: diffusion coefficient and 
reorientation relaxation time 

Fig. |5] shows the density(pressure) dependence of the 
normalized translational diffusion coefficient D/Dq for 
investigated solutes at various temperatures. Here Dq is 
the translational diffusion coefficient of the solute at the 
ambient density for the solvent, which is p = 0.997047 
g/cm 3 for water in cases a), b), and p = 0.782 g/cm 3 
for acetonitrile in the case c), respectively. Fig. |SJi) 
shows that for sufficiently low temperature the diffusion 
coefficient of methanol in water first increases with den- 
sity (pressure), and then begins to decrease. A relatively 



flat maximum is observed for the density slightly smaller 
than that at the ambient condition. The increase in the 
diffusion coefficient with increasing density (or pressure) 
can not be explained by the phenomenological theory 
such as the Stokes-Einstein law. In that reason, such 
behavior is often stated "anomalous" . This is the first to 
realize such anomalous pressure dependence of the diffu- 
sion constant of a solute in water by means of the sta- 
tistical mechanics theory. As can be seen in the figure, 
the anomalous density dependence is suppressed as tem- 
perature increases and turns into normal behavior, i.e., 
the diffusion constant decreases monotonously with den- 
sity. At the ambient condition, the diffusion constant 
does not show any anomaly and is consistent with the 
one obtained from the molecular dynamics simulation 15 , 
which is plotted in the figure with the filled triangles. 
Using the multi-parametric empirical equation of state 
for water— we converted original data 15 from "diffusion 
versus pressure" into "diffusion versus density" . 

Similar features can be observed in the behavior of the 
diffusion coefficient of acetonitrile in water, Fig. |SJd). 
In this case, the increase of the diffusion coefficient at 
the lowest temperature is discernible, but very small as 
can be barely observed with the guide of the cubic spline 
curve. Behavior of the diffusion coefficient at higher tem- 
peratures, as in the case of methanol in water, does not 
show any anomaly. The behavior is in harmony with the 
results of the molecular dynamics simulation at the am- 
bient condition 15 , which is plotted in the figure with the 
filled squares. The experimental measurements of diffu- 
sion coefficients for acetonitrile in water and acetonitrilc- 
D 3 (CD 3 CN) in water, both at T = 303 K, reported 
by Nakahara et al£ testify the same tendency in their 
density (pressure) dependence as those obtained from the 
theory. 

The diffusion coefficient of methanol in acetonitrile are 
plotted in Fig. |5j;) as the function of density. The results 
show monotonous decrease of the diffusion coefficient in 
the entire range of densities and at all investigated tem- 
peratures. 

Fig. shows the density(pressure) dependence of the 
normalized reorientation relaxation time t/tq for inves- 
tigated systems at the same temperatures as for normal- 
ized translational diffusion coefficient D/Dq. tq is the 
reorientation relaxation time of the solute at the ambi- 
ent density for the solvent. Fig. Hjji) is a typical ex- 
ample of anomalous density(pressure) dependence of the 
reorientation relaxation time. It is clear to see that it 
first decreases with density(pressure) and then starts to 
increase. The anomalous behavior becomes weaker as 
temperature increases, and disappears at some point. 
Apparently, the anomaly persists through higher tem- 
perature than what is observed in the diffusion coeffi- 
cient. The theoretical result is consistent again with that 
from the molecular dynamics simulation at the ambient 
condition^ 5 , in the sense that the both plots show a mini- 
mum, although the positions of the minimum are slightly 
different from each other. 
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FIG. 5: Normalized translational diffusion coefficient D/Do 
for methanol in water a); acetonitrile in water b); and for 
methaol in acetonitrile c). □, A, O - theory, ■ A results 
of MD simulation 15 for T = 298 K, lines connecting open 
symbols are cubic splines for the eye-guide. 



FIG. 6: Normalized reorientation relaxation time t/tq for 
methanol in water a); acetonitrile in water b); and for 
methanol in acetonitrile c). □, A, O - theory, ■ A re- 
sults of MD simulation— for T = 298 K, lines connecting 
open symbols are cubic splines for the eye-guide. 



In the case of acetonitrile in water, Fig. IHb), the nor- 
mal behavior of the reorientation relaxation time is ob- 
served almost entirely except for the lowest temperature, 
where a shallow minimum is discernible. The relaxation 
time at the ambient temperature exhibits monotonous in- 
crease. The behavior is in qualitative accord with the re- 
sults from molecular dynamics simulation 15 as indicated 
with the filled squares in the figure. 

Fig. [U) exhibits monotonous increase of the reorien- 
tation relaxation time of methanol in acetonitrile in the 
entire range of densities and at all investigated tempera- 
tures. It can be regarded as a typical example of normal 



dcnsity(pressure) dependence of r. 

Fig. [7|has the purpose to demonstrate the anomalous 
density (pressure) behavior of the reorientation relaxation 
time in terms of the reorientation time-correlation func- 
tion C M (i). It has been found that at T = 298.15 K the 
curves for acetonitrile in water and methanol in acetoni- 
trile move to the right as density increases, namely the re- 
laxation time becomes longer when pressure is increased. 
On the other hand, for methanol in water it shifts first to 
the left, and then for higher densities (pressures) it starts 
to shift to the right. At the highest temperature C M (i) 
for methanol in water moves only to the right, but the 




FIG. 7: Reorientation correlation functions for methanol in water (left), acetonitrile in water (center) and methanol in acetoni- 
trile (right) at T = 298.15 K and different values of density of the solvent. Arrows inside of the each plot indicate direction(s) 
of the displacement for curves with the increase of density. 



displacement is very small. Considering that the time- 
correlation function is related to the reorientation time 
by an integral, the behaviors in the time-correlation func- 
tions are consistent with those of the reorientation time 
as has been described above. 

In our previous studies concerning the neat water, we 
have clarified the cause of anomalous dependence of the 
molecular dynamics on pressure, or the appearance of ex- 
trema in the pressure dependenc o 55 i 58 . The anomalous 
behavior is caused by an interplay of two competing ten- 
dencies: the increasing mobility at the lower pressure and 
the suppression of the motion in the higher pressure. The 
suppression of the motion with increasing pressure is a 
normal behavior caused by increased molecular collisions, 
and the increase in the mobility is caused by the decrease 
in the dielectric relaxation time. The appearance of the 
extrema is attributed to the increasing dielectric motion 
at the lower pressure. The origin of the decrease in the 
dielectric relaxation time at the lower pressure is in the 
reduction of the electrostatic friction on the dielectric 
mode. Although the detailed analysis is presented in our 
previous paper, let us summarize it below. The dielec- 
tric mode describes the collective alignment of the polar 
molecules. When the inhomogeneity of the number den- 
sity is present as is described in Fig. [8^), the alignment of 
the collective dipole moment induces the inhomogeneity 
of the charge density, which leads to the electrostatic fric- 
tion on the dielectric mode. Increasing density or pres- 
sure will suppress the number density fluctuation, which 
in turn causes the reduction of the charge density fluc- 
tuation as is illustrated in Fig. [SJa). The reduced charge 
density fluctuation in turn causes decrease in the elec- 
trostatic friction on the dielectric mode, or increase in 
the collective rotational motion. This is the rough pic- 
ture with respect to the anomalous pressure dependence 




FIG. 8: Schematic view on the origin of electrostatic friction 
on dielectric modes with allowance for large a) and small b) 
density fluctuations. 



of the rotational dynamics in water, which is consistent 
with the analysis based on the mode coupling theory^. 
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Now we turn to the single-particle motion in water. In 
the case of neat water, the translational and the rota- 
tional diffusion of a molecule respond to pressure differ- 
ently, as is observed experimentally. The translational 
motion of a molecule couples with the high wave vector 
collective translational mode, which have inherent ten- 
dency to become slower as pressure increases. On the 
other hand, the coupling between the rotational motion 
of a molecule and high wave vector collective transla- 
tional modes is weak, and the friction on the rotation 
of a molecule mainly comes from the coupling with the 
low wave vector charge density mode, which is usually 
referred to as the dielectric friction. The rotational mo- 
tion of a molecule will induce a relaxation process of sur- 
rounding molecules in order to make themselves aligned 
to the electric field produced by the new orientation of 
the molecule in concern. The energy dissipation asso- 
ciated with the relaxation process is an origin of the 
friction on the rotational motion, or the dielectric fric- 
tion. Slower the dielectric relaxation, greater the dielec- 
tric friction. The anomalous pressure dependence ob- 
served in the translational dynamics or the diffusion co- 
efficient turns out to be the secondary effect induced by 
the anomaly in the rotational motion. The anomaly in 
the translational diffusion is therefore weaker than that 
in the rotational relaxation. 

The solute-solvent systems considered in the present 
paper can be understood in a similar way as follows. 
Whether the anomalous behavior actually shows up or 
not also depends on the role of repulsive core in the ro- 
tational motion of solute, which gives rise to the col- 
lisional friction: a spherical molecule like water will be 
free from the collisional friction, while a rod-like molecule 
such as acetonitrile should have significant collision with 
solvent molecules upon rotation. The behavior exhib- 
ited in Fig. [7| is consistent with the physical picture 
drawn above. Figs. EJi) and0)), which concern water as 
solvent, show anomalous density dependence of the ro- 
tational relaxation time, which disappears with increas- 
ing temperature. The methanol-in-water system exhibits 
more significant anomaly compared to the acetonitrile-in- 
water system. The anomalous behavior is attributed to 
the strong electrostatic interaction, so called "hydrogen 
bond", among the solvent molecules and those between 
solute and solvent. The former causes the decrease in the 
dielectric relaxation time with pressure, and the latter 
brings about the coupling between the dielectric mode of 
the solvent and the rotation of the solute. The anomaly 
is largely suppressed for the acetonitrile-in-water system 
due to the enhanced significance of the repulsive core in 
the molecule. The methanol-in-acetonitrile system does 
not show any indication of the anomalous density depen- 
dence, although there is a hydrogen bond between solute 
and solvent as is indicated in Fig. 0J)). The reason is 
because there is strong collisional friction on the collec- 
tive reorientation of the solvent in this case, so that the 
dielectric relaxation becomes slower with pressure. The 
enhanced density or pressure just amplifies the signifi- 
cance of the repulsive core upon rotation and thereby of 
the collisional friction. 



V. SUMMARY 

In present paper we have calculated the density(pres- 
sure) dependence at various temperatures for the trans- 
lational diffusion coefficients D and rank-1 reorientation 
relaxation times r for three solute-solvent systems at in- 
finite dilution: acetonitrile and methanol in water, and 
methanol in acetonitrile. Calculations have been perfor- 
med using the site-site memory equation with the mo- 
de-coupling approximation for memory kernels, and the 
DRISM theory for static properties. 

The theory was able to reproduce qualitatively all main 
features of temperature and density dependences of D 
and r observed in real and computer experiments. In 
particular, anomalous behavior, i.e. the increase in mo- 
bility with density, is observed for D and r of methanol 
in water at the ambient condition, while acetonitrile in 
water and methanol in acetonitrile show the ordinary be- 
havior, i.e. the monotonous decrease in mobility with 
increasing density. The variety exhibited by the differ- 
ent solute-solvent systems in the density dependence of 
the mobility was interpreted in terms of the two compe- 
ting origins of friction, which interplay with each other as 
density increases: the collisional and dielectric frictions 
which, respectively, increase and decrease with increasing 
density. The cause of the dielectric friction is the inhomo- 
geneity in the charge density of solvent, induced by the 
number density fluctuation due to the hydrogen bond. 
The suppression of the number density fluctuation and 
the consequential decrease in the charge density fluctua- 
tion are the physical origins of the anomalous density de- 
pendence of the molecular mobility in hydrogen-bonding 
systems. 
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